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Abstract 

Here we compute the static potential in scalar QED% at leading order in 1/Nf. 
We show that the addition of a non-minimal coupling of Pauli-type (e^ ua j^d u A a ), 
although it breaks parity, it does not change the analytic structure of the photon 
propagator and consequently the static potential remains logarithmic (confining) at 
large distances. The non-minimal coupling modifies the potential, however, at small 
charge separations giving rise to a repulsive force of short range between opposite 
sign charges, which is relevant for the existence of bound states. This effect is in 
agreement with a previous calculation based on Moller scattering, but differently 
from such calculation we show here that the repulsion appears independently of the 
presence of a tree level Chern-Simons term which rather affects the large distance 
behavior of the potential turning it into constant. 

PACS-No.: 11.15.Bt , 11.15.-q 
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1 Introduction 



An important problem in high energy physics is the lack of a rigorous proof of color 
confinement in 4D QCD. Different techniques have been used to tackle this problem. 
We can mention lattice simulations pQ , supersymmetry j2] and lower dimensional models 

In order to investigate the contribution of the matter fields to this problem we inte- 
grate over such fields in the path integral and derive an effective action for the vector 
bosons, solving the equations of motion of this quantum action we can compute the po- 
tential between two static charges separated by a distance L. A monotonically increasing 
potential as L — > oo signalizes confinement. This route has been followed in [HI El E] in 
the case of QED%. In that model, if we work with two-component fermions, the fermion 
mass term breaks parity and a Chern-Simons term is dynamically generated leading to 
an important change in the analytic structure of the photon propagator which turns the 
classically confining logarithmic potential into a constant at large distances. 

In the case of scalar QED 3 we have a different scenario since its mass term, like the 
rest of the Lagrangian, is parity symmetric and no parity breaking term is dynamically 
generated, so the classical logarithmic potential survives at quantum level. Therefore, it 
is expected that the inclusion of parity breaking terms in the Lagrangian would strongly 
modify the static potential. A natural possibility to be considered is a non-minimal 
coupling of Pauli-type which breaks parity but preserves gauge invariance. This term is 
rather simple in D = 3 where the dual field strength (F^ = e llva A ii d v A a ) is a pseudo- vector 
and the addition of the non-minimal coupling amounts to the replacement eA^ — > eA^ + 
7-F M where 7 is the non-minimal coupling constant which has negative mass dimension. 
This term has been considered before in the literature of QED 3 and scalar QED 3 , see 
e.g. [9-16]. 

Another motivation for the inclusion of the non-minimal coupling comes from ^21 
[T3J E] where there are indications, see however jT3] , that the coupling of a gauge field 
to fermions via a Pauli term could give rise to anyons with no need of a Chern-Simons 
term. Since the change of statistics is a long range phenomena and the Chern-Simons 
term indeed changes the static potential at large distances, we would like to include the 
Pauli-type interaction in order to check, at least in some approximation, if it could really 
produce large distance effects. 

A further point concerns previous calculations in the literature. It has been claimed 
in jini E| that the effect of the non-minimal coupling on the static potential only appears 
if a Chern-Simons term is present. This is apparently not the case of QED 3 with four- 
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component fermions where no Chern-Simons term is generated but still there is some 
influence of the Pauli-type term on the static potential at low distances [9j. It is important 
to remark, however that here and in |5| one works at leading order in 1/Nf which requires 
the calculation of the one loop vacuum polarization diagram, while the calculations of 
|lbl IT7] are based on the one photon exchange diagram at tree level (Moller scattering) 
in the non-relativistic limit. In order to control the effect of the Chern-Simons term 
and compare our results to [H)J E] we introduce here, besides the Pauli-type term, a 
Chern-Simons term at tree level with an arbitrary coefficient. 

We have already mentioned that the Pauli term demands a coupling constant with 
negative mass dimension (non-renormalizable) so we found suitable to use 1/Nf expan- 
sion since there are some arguments ^H] m favor of the 1/Nf renormalizability of such 
interaction. In the next section we start by integrating over the Nf scalar fields at leading 
order in 1/Nf. Then, we analyze the analyticity properties of the corresponding photon 
propagator. In section III we minimize the effective action and compute the static poten- 
tial V(L) numerically for a finite scalar mass and analytically in the limit m — > oo. We 
draw some conclusions in section IV. 

2 The photon propagator at Nf —> oo 

Our starting point is to integrate over the Nf scalar fields <p r , r = 1, 2, ■ • • , Nf in the 
partition function below: 

Z = fvA,e l f d " x J] V € V^r e~i 1 ^ W]</v 

= C JvApj! rf3 *£(^icx t )-^/Trln[DM DM+m 2] ^ 

where C is a numerical constant and 

£ (^>iext) = -~ F% ~ ~e^d»A a + | {d^f - A u3 ^ t (2) 

The external current corresponding to a static charge Q at the point (xi,X2) = (L/2,0) 
is given by j" xt = Q8(x2)S(xi — ~)8 v0 . Later on, the interaction energy of a cou- 
ple of charges — Q and Q separated by a distance L will be calculated via V(L) = 
—QAq (x\ = —L/2, X2 = 0) where Aq(x u ) will be obtained minimizing the effective action 
coming from The covariant derivative: D^cp = d^4> — iecfyA^/ \fN~f — 27 4>F^/ \jNf 
includes the non-minimal coupling constant 7 which has negative mass dimension [7] = 
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—1/2 while [e] = 1/2 and [9] = 1. The dual of the strength tensor is defined here as 

rp , f)v A a 

r fi tfiuaU Si ■ 

The next step is to evaluate the trace of the logarithm perturbatively in 1/Nf. We have 
two types of interaction vertices coming from C^ t = % (0*<9 M — 0<9 M 0*) [eA^ + ^F^) / \/Nf 

( 2) 2 

and C\J t = i(f)*(f) (eA^ + ^F^) /Nf. Thus, the leading contribution in l/Nf would come 
from just one vertex of the first type, however, since it involves derivatives of the scalar 
fields the Feynman rules for scalar QED include a factor p™ + p° ut where those are the 
incoming and out-coming momenta of the scalar fields. Therefore the diagram (tadpole) 
will be proportional to the integral J <i 3 p (p™+p° ut )/(p 2 — m 2 ) = J d 3 p2p tJ- /(p 2 — m 2 ) which 
vanishes in the dimensional regularization adopted here. The next leading contribution 
includes either two vertices of the first type or one vertex of the second type. Both con- 
tributions will be independent of Nf due to the overall factor Nf in front of the logarithm 
in (0) and will survive the limit Nf — > oo. The next contribution would come from three 
vertices of the first type and would be of order 1/ v/iV/ so it vanishes if Nf — > oo. Such 
contributions and higher ones will be neglected henceforth. In conclusion we have, up to 
an overall constant, Z = J T>A^e tScS where: 

S eS = J d 3 xC - % - J Sk (eA»(k) + iHk)) T, u (eA»{-k) + ^{-k)) (3) 

The quantities A^, F v are Fourier transforms and 

T Q/3 = -2^/ (1) + ig (4) 
Using dimensional regularization we have obtained for the Feynman integrals: 

w= f d 3 p 1 m 

J (2vr) 3 p 2 - w? 4tt 1 ; 



T (2)_ f d 3 p (2p + k) a (2p + k) /3 vm _ 0yf n i 

^ ~ J (2vr) 3 (p 2 - m 2 )[(p + k) 2 - m 2 ] ~ 8vr [ 9aP H ^ 

With z = k 2 /4m 2 and 9 af3 = g a/3 — k a kg/k 2 . In the region < z < 1 we have 
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Above the pair creation threshold (z > 1) the integral 1^1 develops a real part which 
will be neglected henceforth. For future use we have given the large mass expansion for 
f 2 . The static potential requires the expression for the effective action for z < which 
can be obtained by analytically continuing ((7j) and (jHI). Namely, with z = —z > we 
have : 



h 
A 



-fi 



—= arctan v z 



z H z 

3 15 35 



(9) 
(10) 



Our result for T a p is in agreement with JH] and it is transverse k a T a p = = T a pk^ in 
accordance with gauge invariance. Now we can write down the effective action for scalar 
QED 3 including vacuum polarization efffects : 



S, 



cff 



d 3 x 



.1 

4 



7 2 D/ 2 , e 2 f 2 



+ 



lQirm lQirirr 



-e 



fiua 



A»d u A a - 



e7 



16m7r 



€„ va A»d v Uf 2 A a - A v j, 



V 

cxt 



where f 2 = D/4m 2 ) is given in (JZj) and ©• Notice that, besides the tree level Chern- 
Simons term, another parity breaking term appears in (jllj) due to the magnetic moment 
interaction. Although the action (fTTj) is non-local it can be made local in the large mass 
limit m — > oo as in [20J. Introducing the dimensionless constants 



e 2 _ ej 9 

C\ — 5 c 2 = — 5 C3 — 

167rm 87r ' 2m 



(12) 



Taking m — > 00 while keeping the dimensionless constants finite, the only effect of the 
vacuum polarization is a finite renormalization of the Maxwell term, i.e., 



S P fe (m — > 00) 



ft 



(13) 



On the other hand, for finite mass we can write down: 



S, 



cff 



d 3 xd 3 y 



A»{x) D ^ V) A»{y) - A u f cxt 5^(x - y) 



rf3fc - A"(-Jfc)- / ft^jl 



(2vr 



(14) 
(15) 
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Where the photon propagator in momentum space is given by 



D 



a (9, 



.0 



(16) 



with 



1 



(17) 




b 



Cl (D + + D_) 



(18) 



8m 2 y/zD + D 



.ci (D+-D_) 



(19) 



c 



— z 



16m 3 zD,D_ 



D± = \/z ci + (ci ± \fzc 2 ) f 2 ± C1C3 = g±(z) ± C1C3 



(20) 



Now we are able to analyze the analyticity properties of the photon propagator. First of 
all, we notice that the massless pole z = in the denominator of which is typical of 
a Maxwell- Cher n- Simons theory, is a gauge artefact. It disappears from gauge invariant 
correlators involving the field strength F^. It can be shown [21J to have a vanishing residue 
(non-propagating mode). Since the factor \J~z in the denominator of (JTSJl is cancelled ed 
out by the numerator, the only possibilities for poles in the propagator stem either from 
D + = or D_ = 0. Due to c\ > and 2/3 < f 2 < 1 we have g±(z) > and consequently 
we can only have D + = or = for C3 < or C3 > respectively. We never have 
two poles at the same time. In the absence of the Chern-Simons term, i.e., C3 = 0, the 
product D + D_ will be proportional to z and we are left with one massless pole z = 0. 
Since \im z ^o z(a + b) = —1/ (1 + 2ci/3) < the residue at this pole will be positive 
and this represents a physical massless photon which will be responsible for a long range 
logarithmic static potential. On the other hand, if C3 7^ 0, since the denominator D + D_ 
is symmetric under c 2 — > — c 2 ; c 3 — ► — c 3 , it is enough to consider only D + = assuming 
C3 < the other case D- = with C3 > follows from the symmetry. Numerically, 
we have checked that whatever sign we choose for c 2 the function g+(z) is monotonically 
increasing and satisfies g+(z) > consequently its maximum is <7+(l). Therefore, see (|2*0j) . 
if c 3 < — g + (l)/ci = [(ci + c 2 ) 2 + Ci] /ci then we have no poles and so no particle in the 
spectrum. On the opposite, if —g + {l)/c\ < C3 < we are always able to find numerically 
one massive pole for some < z < 1 such that D + = which is a typical effect of a 
Chern-Simons term [22j- As we move toward the left limit value C3 — ► — g + (l)/ci, the 
photon mass increases to the point where it reaches the pair creation threshold k 2 = 4m 2 
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at c 3 = —g + (l)/cx. Due to the symmetry c 2 — > — c 2 ; c 3 — > — c 3 we conclude that whenever 
the tree level Chern-Simons term is present and its coefficient is not too much negative 
or too much positive ( | C3 1 < gr + (l)/ci) we have one massive physical (positive residue) 
photon and if the Chern-Simons term is absent we have one physical massless photon. 

It is remarkable to find a no poles region in the photon propagator. One might think 
that this is due to some convergence problem of the 1/Nf expansion which has been 
introduced because of the non-renormalizable non-minimal coupling c 2 . However, even if 
c 2 = the Chern-Simons coefficient must obey an upper bound |c 3 | < g+(l)/ci = 1 +Ci in 
order to have a physical pole in the photon propagator at one loop level. By analytically 
continuing, see (J5j). the expression for the propagator to the region z = k 2 /4m 2 < 
we have checked that tachyons can only appear for a special fine tuning of the coupling 
constants for which we did not find any special interpretation, namely, the tachyonic 
pole must be a solution of zf% = c 3 /(2c2) and this solution must be such that cfc 3 = 
— 5c 2 (2ci — c 2 c 3 ), although explicit numerical solutions are possible we have found those 
fine tuned cases rather artificial. In particular, they have apparently no relationship with 
the no-pole region (|c 3 | > g + (l)/ci) and will be disregarded in this work. In the next 
section we use the photon propagator as an input to calculate the static potencial V(L). 

3 The static potential V(L) 

Minimizing the effective action ()14|) we obtain: 



Mv) = j d 3 xDp a (y,x)j: xt (x) (21) 
= / ^^(k) j d 3 xe^y^ 3 : xt (x) (22) 

Since the external current is time independent, in (|22|) there will be a factor J dxoe~ lk ° x ° = 
2%5(ko) which allows an exact integration over fco, implying k^k^ = —k\ — k\ = —k 2 < 0, 
consequently z < and we need the analytic continued functions / 2 instead of / 2 . The 
integrals over x\ and x 2 can also be readily done using the delta functions in the external 
current. The angle part of the integral dk\dk2 = kdkdO gives rise to the Bessel function 
J (kL). Thus, we are left with the radial integral over k = \/kj + k\. Placing the negative 
charge — Q at (yi,y 2 ) = ( — L/2,0) we have 



V(L) = -QA ( yi = -L/2,y 2 = 0) 
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O 2 f°° 

= -— / dkkbJo(kL) (23) 
27r Jo 

The tilde in the expression b stands for the analytic continuation of (JTHJl to z < 0. 

Now we discuss some special cases starting with the pure scalar QED3 where C3 = = 
c 2 . In this case kb J (kL) = J (kL)/ k{l + C1/2) since f 2 (k = 0) = 2/3 and Jq(0) = 1 
we have an infrared divergence at k = and the integral (|2*3*|) is divergent as it stands. 
We make a subtraction in order to get rid of this infrared divergence and define: 

V(L) - V(L ) = - — lim / dkkb [J (kL) - J (kL )} (24) 

27T x-*0 J x 



In general the integral (|24|) must be calculated numerically, one exception is the large 
mass limit m — ► 00. In this case kb — > 1/ [k(l + 2ci/3)] and the integral can be calculated 
exactly [2*!?] : 

where 



[n£)-n£o)] m ^ = ^ln^-] (25) 



L ~ t ~ 247i-mJ 

The classical potential is given by (|23jl with Qr replaced by the bare charge Q. There- 
fore, the sole effect of the vacuum polarization is a finite renormalization of the charge. 
The situation is similar to QED3 with four-component fermions where no Chern-Simons 
term is dynamically generated, the only difference is the renormalized charge which is 
Q/ [1 + e 2 /(67rm)] instead of (J2EJ)- Thus, the renormalization factor is larger for fermions 
than for scalars. For finite mass the potential must be calculated numerically. We plot 1 
the results in figure 1 for the masses m = 1 and m = 0.01 and compare with the clas- 
sical result and the result of [H] for four-component fermions. We notice that the finite 
renormalization due to the vacuum polarization is always stronger for fermions than for 
scalars and its effect increases with the mass of the matter fields. For both scalar QED% 
and QED% we see in figure 1 that the numerically calculated static potential at m = 1 is 
already very close to the analytic result (solid lines) obtained in the limit m — ► 00. 

Next, we check the effect of the non-minimal coupling c 2 7^ in the absence of the 
Chern-Simons term (03 = 0). In figure 2 we see that for L —>■ 00 the effect of the 
non-minimal coupling in the vacuum polarization disappears and the potential becomes 
1 In all figures in this work the symbol V stands actually for the difference V(L) — V(L ) 
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Figure 1: The static potential for pure scalar QED% (dark dots) and pure QED% with four- 
component fermions (light dots). The dashed line corresponds to the classical potential. 
We have fixed (ci, c 2 , c 3 ) = (1, 0, 0) and m = 0.01 for the two dotted curves closer to the 
classical potential while we have m = 1 for the farther curves which overlap with the 
m — ► oo analytic result (solid curves) 

logarithmic which can be explained technically by the fact that the Bessel function J (kL) 
oscillates with decreasing amplitude as L — > oo and so the integral will be dominated by 
the pole at the origin k = which makes the higher derivative, see ([lljh contribution 
of the non-minimal coupling negligible. However, in a finite range close to L = the 
non-minimal coupling gives rise to a surprising repulsive force in a much similar way to 
what happens in the case of four-component fermions in [§]. Such repulsive force may 
play an important role in the existence of bound states. Differently from the calculation 
based on the Moller scattering [IB] we see here effects of the non-minimal coupling even 
in the absence of the Chern-Simons term. 

Now we turn on the Chern-Simons term (c 3 ^ 0). As we see in figure 3, the potential 
V(L) — V(Lq) tends to the constant — V(L ) as L — > oo like the case of QED 3 with two- 
component fermions, see 6J, where a Chern-Simons term is dynamically generated. Once 
again, a repulsive force appears for small separations as we switch on the non-minimal 
coupling. As L —>■ oo the only effect of the non-minimal coupling is to change the constant 
—V(Lq). Although, the plot in figure 3 has been made for m = 1 and C\ — 1 we have 
checked that the same form of the potential persists for other values of those constants. In 
summary, the effect of the non-minimal coupling is qualitatively the same in the presence 
of a Chern-Simons term. 
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Figure 2: The static potential for scalar QED 3 without the tree level Chern-Simons term 
(C3 = 0). The solid line corresponds to pure scalar QED 3 (c 2 = 0) while the dark (light) 
dots to c 2 = 2(c 2 = 4). We have assumed m = 3 and c x = 1. 



v 




-i 

-1.5- 



Figure 3: The static potential for scalar QED 3 with the tree level Chern-Simons term 
(c 3 = 1). From the lightest to the darkest curve we have c 2 = 0; 0.4; 0.8. We have assumed 
m — 1 = c\. 
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4 Discussions and Conclusion 



In the case of QED% with two-component fermions it is well known that a Chern-Simons 
term is dynamically generated which makes the photon massive and turns the classical 
confining logarithmic potential into constant at large distances. In scalar QED 3 the 
mass term for the scalars is parity invariant and no Chern-Simons term is dynamically 
generated and so the classical logarithmic potential survives vacuum polarization effects. 
Here we have explicitly confirmed that fact and analyzed the effect of adding a parity 
breaking non-minimal coupling term of Pauli-type as well as a tree level Chern-Simons 
term. It turns out that the non-minimal coupling by itself neither affects the analytic 
properties of the photon propagator nor changes the large distance behavior of the static 
potential which is unexpected from the point of view of the interpretation that this term 
may originate anyons with no need of a Chern-Simons term, see |12[ IT4l IT3] but see also 
|lfij . On the other hand, at small charge separations the non-minimal coupling gives 
rise to a repulsive force between opposite sign charges which has been observed before 
in by computing the one-photon exchange diagram (Moller scattering) and taking 
the non-relativistic limit. Notwithstanding, the effect found in |TBJ only appears in the 
presence of the tree level Chern-Simons term and its attractive or repulsive nature depends 
on the sign of (1 — j9/e) contrary to what we have found here where the non-minimal 
coupling influence is present, see figure 2, even if 6 = and its effect is always repulsive 
independently of the sign of 7 or 6. 

Concerning the tree level Chern-Simons term, as expected, it gives mass to the photon 
and shifts the zero momentum pole in the integral involved in the static potential (j24)l . The 
absence of a singularity in the integration path allows us to take the limit L — > 00 before 
performing the integral and so it will vanish as a consequence of Jq(x — ► 00) — > 0. This 
effect of the Chern-Simons term was certainly not surprising. However, it is remarkable 
that we found an upper bound for the absolute value of the Chern-Simons coefficient in 
order to have a physical pole in the photon propagator at one loop level. As we increase 
such absolute value the photon mass increases and penetrates the real pair creation region 
k 2 > 4m 2 for finite values of the coupling constants of the theory. We can mention that 
this situation is not peculiar to scalar fields since we have noticed in |H] that it happens also 
in QED3 with two-component fermions. In that case, if c\ = e 2 / (Wirm) > 1 there will be 
no poles in the photon propagator at one loop level. However, one could argue that c\ is a 
dimensionless constant which controls the perturbative expansion (for Nf = 1) such that 
the upper bound could be understood as a limit for perturbation theory. This argument 
does not work for scalar QED 3 even if we drop the non-minimal coupling (02 = 0) since 
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the upper bound increases with c% which makes the latter case more intriguing. 

At last, we notice that the static potential in pure scalar QED% without tree level 
Chern-Simons term has been studied in [21] where the authors conclude that the potential 
is of screening type and even fractional charges can be fully screened. However, the 
authors of [21] have neglected terms of order e 2 /m which have been considered here. 
Besides, they have gone above the pair creation threshold and made use of variational 
methods altogether with a peculiar Ansatz for the two particle wave function, so we can 
hardly compare their findings with our results obtained below the pair creation threshold. 
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